A close look into the excluded volume effects within a double layer 
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We explore the effect of steric interaction on the ionic density distribution near a charged hard wall. 
For weakly charged walls, small particles, and monovalent ions the mean-field Poisson-Boltzmann 
equation provides an excellent description of the density profiles. For large ions and large surface 
charges, however, deviations appear. To explore these, we use the density functional theory. We 
find that local density functionals are not able to account for steric interactions near a wall. Based 
on the weighted density approximation we derive a simple analytical expression for the contact 
electrostatic potential which allows us to analytically calculate the differential capacitance of the 
double layer. 

PACS numbers: 



I. INTRODUCTION 

The standard Poisson-Boltzmann equation (PB) is 
accurate for dilute aqueous solutions containing small 
monovalent ions. The first approximation of the PB the- 
ory is the mean-field account of the electrostatic interac- 
tions, where the discrete nature of ions is neglected and 
each ion interacts with the mean-potential generated by 
the local charge density distribution p c (r) according to 
the Coulomb law, ^i(r) = (4-Tre)- 1 / dr' p c (r')/\r' - r|, 
where e is the dielectric constant. This description 
clearly ignores ionic correlations. The other approxima- 
tion within the PB theory is the lack of internal structure 
of each ion: an ion is a point in space characterized only 
by its charge. Under this reduced description the ions 
F~, Cl~, and NOJ" are indistinguishable, although their 
size [l[ and polarizability are obviously different. The 
PB description, therefore, breaks down when the above 
reductions can no longer be justified. In the present work 
we consider the weakly correlated limit so that the mean- 
field description of electrostatics holds, but we go beyond 
the standard PB equation by including the nonelcctro- 
static excluded volume effects [l|. In aqueous solutions 
the excluded volume interactions are further increased on 
account of hydration, the binding of water molecules to 
ions, so that the effective ionic radii are larger than the 
crystallographic ones . Another effect of hydration, 
related directly to electrostatics, is the dielectric decre- 
ment [fUl] associated with the decrease of the medium di- 
electric constant on account of polarization saturation of 
a solvent constituting the hydration shell. In the present 
work we do not consider this effect and focus strictly on 
the non-electrostatic excluded volume effects. 

The excluded volume effects reorganize the double 
layer so that a density profile is no longer monotonically 
decreasing, but it acquires oscillatory structure reflect- 
ing molecular composition of an electrolyte @. Clearly, 
the modifications of the density profile lead to modified 
electrostatic properties. In the present study we take a 
close look into the excluded volume effects where ions 
are represented as hard spheres with the same diameter, 



the so-called, restricted primitive model. Not only do 
we investigate how the ionic system responds to these 
effects under various conditions, but also how well var- 
ious theories reproduce this behavior. The hard-sphere 
interactions lead to the nonlocal effects, demanding non- 
local functionals to account for the oscillatory structure 
in the density profile. The nonlocal treatments, such as 
the density functional theory (DFT), can be computa- 
tionally demanding if the system is sufficiently complex. 
Frequently in electrostatics, a variant of the local den- 
sity approximation (LDA) is used, where the ideal gas 
entropy is substituted by that of the lattice-gas lead- 
ing to the modified Poisson-Boltzmann equation (MPB) 
dHHH. With growing popularity of the MPB as a simple 
and tractable model, we feel that a careful and system- 
atic study of this method is still lacking. First of all, 
the MPB is a local theory inappropriate for moderately 
and highly inhomogeneous fluids. Furthermore, even as a 
theory of a homogeneous fluid it already fails at predict- 
ing the correct second virial term. We consider another 
variant of the LDA based on the scaled particle theory 
(SPT). Surprisingly, although highly accurate for treat- 
ing homogeneous fluids, it is less accurate than the MPB, 
which is unreliable for homogeneous fluids. We compare 
all the results against the nonlocal fundamental measure 
density functional theory [l2j|. Our main focus is the re- 
gion near and at the contact with the wall, where exact 
thermodynamic relations are known. We investigate the 
scaling properties of the contact quantities and propose 
an analytical fit to the DFT data that permits an accu- 
rate prediction of these quantities, without resorting to 
numerical computations involving the entire region. 

Throughout the paper we make frequent use of the ini- 
tials: DFT, LDA, MPB, and SPT. To avoid confusion we 
give a quick account of each. DFT always refers to the 
fundamental measure density functional theory as origi- 
nally derived by Rosenfeld [12(- This is the only nonlo- 
cal theory used in this work. LDA refers to the general 
class of approximations where a density functional is de- 
fined locally based on a free energy of a homogeneous 
fluid. MPB is the special example of the LDA, when the 
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excluded volume effects are represented as that of a ho- 
mogeneous lattice gas [HI- And the SPT stands for the 
scaled particle theory and refers to the theoretical ap- 
proach based on scaling a particle size fl3l - [l6l ]. We use 
SPT to refer to another LDA approximation based on 
the SPT results for a uniform fluid. In the paper two 
LDA approximations will be considered: the MPB and 
the SPT. 

The paper is organized as follows. In Section II we dis- 
cuss a reversible work of inserting a test particle into a 
hard sphere fluid, W^ a , and emphasize its nonlocal char- 
acter by deriving the exact result for inserting a point 
particle into a hard-sphere fluid. We then introduce 
the fundamental density functional theory as an accu- 
rate nonlocal treatment of the hard-sphere interactions. 
In Section III we review the contact value theorem and 
explore its implications for the counterion distribution 
near a hard charged surface. We point out that the LDA 
theories lead to contact relations different from the exact 
contact value theorem. This is the result of represent- 
ing the overcrowding in the double-layer as the density 
" saturation" . In Section IV we investigate the scaling of 
the quantity A w , which represents the reversible work of 
brining the uncharged fluid particle from the bulk into 
contact with the wall. We propose the analytical for- 
mula for A w that fits the DFT data points. In Section V 
we apply this formula to calculate the surface potential 
and the differential capacitance. Section VI finalizes this 
work with the concluding remarks. 

II. PRELIMINARIES 



with respect to a density distribution, ^ = 0, and leads 
to 

Pj (r) = pge-^M+^M-^l. (5) 

where ip = |^ /qj is the mean electrostatic potential, 
Whs = g^i and p b j is the bulk density of an ion species 
3- 

Physically, Wba represents the reversible work of in- 
sertinga test particle at position r into a hard-sphere 
fluid j!3l - [lq ] . If the system has a wall, the insertion per- 
formed far from the wall, Wh s (°o) = W^ s is identical 
with the excess (over the ideal gas) chemical potential. 
Alternatively, Wh s can be seen as a depletion interaction: 
as a particle approaches a hard-wall, it feels itself being 
pushed towards the wall by the other particles. Appli- 
cation of the Poisson equation, eV 2 ip = —p c , leads to a 
kind of modified Poisson-Boltzmann equation, 

eVV = -J2 p)q j e-^^ +w ^- w ^\ (6) 
3=1 

where the surface charge a c is determined by the bound- 
ary conditions at a charged wall, 

dip 
ox 

In order to obtain the potential, one still needs some 
kind of closure for Wh s that will determine the excluded 
volume effects. 



We write the free energy functional as having three 
separate contributions, 

F = Fid + F c + Fhs, (1) 

the ideal gas contribution, 

K 

F id = k B Tj2 / drp,(r)(log ft (r)A 3 - 1), (2) 

3=1 

the mean-field electrostatic interactions, 

F c = — f dr [ dr > [Pc(r)+p f (r)][p c (rQ+p f (r')] 
8ire J J |r' - r| 

(3) 

and the excluded volume contributions, F^s, where all 
ions are taken to have the same diameter a. pi denotes a 
fixed charge of a macromolecule and in the present work 
we consider a charged wall, Pf = a c 5(x), where a c is the 
surface charge. p c = Ylj=i I3P3 ' IS the charge density of 
mobile ions, where qj is the charge of an ion species j and 
K is the number of all ionic species. Finally, A is the de 
Broglie wavelength. The equilibrium density minimizes 
the grand potential functional, 




A. insertion of a point particle — the exact case 

IVhs is a nonlocal quantity. To see the origins of the 
non-locality we consider a simple case where the insertion 
work can be expressed as an exact density functional. 
The case we refer to is the reversible work of inserting 
a point test particle into a hard-sphere fluid, W p . The 
range of the hard-core interaction between the test point 
particle and the fluid hard-spheres is half the diameter, 
a/2. The reversible work is related to the insertion prob- 
ability, lip = e~ m v. For a given instantaneous configu- 
ration of a hard sphere fluid, the space is either occupied 
or is empty. The point particle can be inserted into the 
empty space between the hard spheres. The instanta- 
neous insertion probability, n p , is either zero, if we hit 
on an occupied space, or one, if a cavity is found, 

jv 

n p (r) = l-^0(|-|ri-r|) 

1=1 

= 1- J dv'p{T')6(^-\v'-v\), (7) 

where p(r) = Yli=i ^( r » — r ) ^ s the density operator, N is 
the total number of fluid particles, and 0{x) is the Hcavi- 
sidc step function. Averaging over all the configurations, 
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we obtain 

(n p (r))=n p (r) = l-| dr'p(r')9(~-\r>-r\). ( 8 ) 

The reversible work of inserting a point particle at posi- 
tion r is 



pW p (r) = -log 



1- Jdr'p{T')0{^ 



•I) 



(9) 



The non-locality of W v is already evident in this simple 
example and shows that the work of inserting a point par- 
ticle depends on the weighted average of the inhomoge- 
neous density distribution. The extent of the nonlocality 
is determined by the range of the hard-core interaction. 
To obtain an expression for Whs, & reversible work as- 
sociated with expanding a point particle to diameter a 
has to be considered, and W p constitutes only a part 
of the entire quantity Whs- For the homogeneous fluid, 
this procedure leads to the scaled particle theory (SPT) 
[l3j . in which the work of expansion to finite diameter 
a is interpolated between the two exact limits, A = oo 
and A = a/2, where A is the range of the hard sphere 
interaction between a test and a fluid particle. For an 
inhomogeneous fluid the derivation is more involved fl2| . 



B. density functional theory 

The ideas of the SPT applied to inhomogeneous flu- 
ids lead to the formulation of the fundamental mea- 
sure density functional theory (DFT) [jj]. The DFT 
approximates the excess free energy due to the hard- 
sphere interactions through the free energy contribution, 
/3F hs [p\ = J dr$ hs , where 



$hs = -polog(l-/0 3 ) 



PlP2 - Pi ■ P 2 , P2 



Pi 



P2 ■ P2 



1 



P3 



24ir (l-p 3 ) 



(10) 

The hard sphere contributions are locally represented 
through weighted densities, instead of a physical density, 
as is the case for the LDA approximation, 



J dr'p(r')u a (r-r') (11) 



with w Q (r — r') denoting a weight function [12],|l7j. There 
are five weight functions which are obtained through de- 
composition of the Mayer-/ function into convolution 
products [HI . The two of the weight functions are vector 
quantities and represent a distribution of normal to the 
sphere surface vectors. In the uniform fluid they have no 
contribution. The weight functions arc given below, 



u)2v(r) = ~5{R-r), 
r 



wa(r) = S(R-r), ws(r) = 9(R-r), 

r 

(12) 

and Wo(r) = w 2 (r)/(47ri? 2 ), U\(t) = o; 2 (r)/(47ri?), and 
u>i V {v) = W2v( r )/(47r_R), where R = a/2. Since we as- 
sume the size of all ions to be the same, the weight func- 
tions are the same for all the ionic species. The insertion 



work Whs is defined as the functional derivative of the 
excess part of free energy, 



PWhs{r) 



S(3F hs 
Sp 



E 



, d$ hs (r') , 

dr — — uj a {r -r). (13) 

op a 



The DFT theory presented above accurately captures the 
dimensional crossover when the 3d fluid is confined to a 
quasi- 2d space [l9| — if a degree of freedom is completely 
eliminated in the transverse direction, a 2d fluid is accu- 
rately recovered. 




FIG. 1. The distribution of counterions near a wall. The 
points are the results of the Monte Carlo simulation of 1 : 1 
electrolyte between two parallel charged walls with the surface 
charge a c = 0.1C/m 2 and the separation L = 3.52nm. The 
remaining parameters are: the particle diameter a = 0.4nm, 
the salt concentration c s = 9.5M, and the Bjerrum length 
As = e 2 /eaksT = lnm, where e is the elementary charge. 
The solid line is the DFT prediction. For technical details 
regarding the simulation see the reference [2C| . 



C. local density approximation (LDA) 

Although highly accurate, see Fig. ([1]), the numer- 
ical implementation of the fundamental measure den- 
sity functional theory is quite involved. To simplify the 
calculations, local density approximations (LDA) have 
often been invoked. Within the LDA, the thermody- 
namic relations derived for a homogeneous fluid are ap- 
plied locally to an inhomogeneous situation. For ex- 
ample, within the LDA the insertion probability Whs 
is supposed to be a function of the local density p(r) 
and is identified with the "local" excess chemical poten- 
tial, Whs(r) = Whs(p(r)). Clearly the specific form of 
the LDA will depend on the approximate form of Wi ls - 
The MPB [ll| relies on Wh s derived from the lattice-gas 
description of electrolyte, pWhs = — log(l — 77), where 
r/ = ira 3 p/6 is the packing fraction. The SPT expression 

is = iog(i - n) + 14 " 2 -(i% a )t 5 " 3 - 
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III. THE CONTACT VALUE THEOREM 

In this work we focus on the region at a contact with 
a charged wall. In this region the excluded volume ef- 
fects are most pronounced since the ion concentration is 
highest, and for many problems the contact quantities 
are of primary interest. This is the case for calculation of 
the differential capacitance [2l| , or the estimation of the 
strength of the electrostatic interactions between colloids 
trapped at an interface [22j ■ The theoretical investigation 
in this region is facilitated by the existence of exact ther- 
modynamic relations. This section reviews some analyti- 
cal results valid at a contact with a wall. Throughout the 
paper we make frequent use of either the subscript or the 
superscript w to indicate that a quantity is taken at con- 
tact with the wall. In the same way we use b to indicates 
bulk values. We remind that the charge and number den- 
sity are p c = Y,f=i VjPj and P = Jlf=i Pit respectively. 
Consequently, the bulk and contact value of the number 



density is pb = J2f=i P b j an( l P 



2^,=1 Pj 



respectively. 



exact results 



The simplest expression of the contact value theorem 
relates the bulk pressure P to the average momentum 
transferred at the hard-wall: k B T p w = P, where p w is 
the fluid density at the contact with the wall. If the hard- 
wall has a surface charge a c and the particles are hard- 
spheres with a central charge, then the contact value the- 
orem takes the form [2J , 



k B Tp u 



P+^. 
2e 



(14) 



Within the DFT formalism P = k B Tp b ( 1 + rj b + 1}\ ) / ( 1 - 
rj b ) 3 the pressure is given according to the scaled particle 
theory. This result is a highly accurate approximation 
and the contact value theorem according to the DFT is 
essentially exact. The contact value relation is the result 
of the mechanical force balance condition, 



representing the reversible work expanded to bringing a 
hard-sphere particle (with no electric charge) from the 
bulk electrolyte to the wall. This is not to say that A w 
is independent of electrostatic properties. The work of 
bringing an uncharged hard-sphere to the wall depends 
on the ionic distribution within the double-layer that de- 
pends on the electrostatic interactions. Substituting the 
density in Eq. (fTB"]) into the contact value theorem gives 



k B T 



K 

E 



= P 



2e 



(18) 



For the case A w = 0, ip w is known exactly and is de- 
termined by the bulk quantities and the surface charge. 
The standard PB equation describes this scenario. A fi- 
nite A w modifies the contact potential. A w > signifies 
that it is easier to insert a particle far away from the wall 
than at the contact. This is intuitive because near the 
wall counterions are in excess. Consequently the contact 
potential is larger than that obtained from the standard 
PB equation — steric interactions deplete some ions from 
the first layer making it more difficult to neutralize the 
wall charge. On the other hand, A w < signifies that it 
is easier to insert a particle at the wall than far from it. 
Although this seems counterintuitive, it can occur at low 
values of the surface charge a c and/or large rj b . To see 
this we consider the limit a r — > 0, 



lim Ai 



log 



p 



k B Tp b 



(19) 



In this limit, A w is always negative (since P > k B Tpb), 
and the larger the bulk packing fraction rjb, the larger 
will be the modulus of A w . At large ionic concentrations 
- such as, for example, found in ionic liquids - the range 
of <t c where A w is negative can be significant. In Fig. ([2]) 
we plot A w as a function of rjb for different values of the 
surface charge. The crossover where A„, changes sign 
to negative is pushed to larger rjb as a c becomes larger. 
Eventually the crossover is suppressed as rjb passes the 
freezing transition (at rjb ~ 0.49). 



pV [k B T log p] = p c V4> + P VW hs , 



(15) 



where the thermodynamic force on the left-hand side is 
counterbalanced by the electrostatic force and the deple- 
tion force due to the collisions with other fluid particles 
on the right-hand side. 

The explicit expression for the contact density, using 
Eq. ©, is 



K 
3=1 



b e -Pqjip u 



-A, 



where we have introduced the quantity 



A, 



dx 



dWhs 
dx 



(16) 



(17) 



B. local density approximation (LDA) 

Within the local density approximation the insertion 
work is a local quantity, 



IF, 



9$ 



hs 



lis 



dp 



P=p(r) 



(20) 



and the force density due to the hard-sphere interactions 
can be related to the gradient of a local pressure, 



- P VW hs = - V 



P^\ *hs 

dp 



= -V[P-k B T Pb }. (21) 



The force balance equation in Eq. (jT5j) becomes 

Pc ViP = -VP, (22) 
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FIG. 2. The reversible work of dragging an uncharged test 
particle from bulk to the wall plotted as a function of the 
bulk packing fraction. The data points are from the nu- 
merical DFT calculations for a symmetric salt 1:1. The 
bulk ion concentration is pb = 2c s , the bulk packing fraction 
77b = 7rcr 3 c s /3, where c s is the bulk salt concentration. The 
parameters are: the particle diameter is a — 0.8nm, the Bjer- 
rum length is As = 0.8nm, and the bulk salt concentration 
is c s = 0.1M. Negative A w signifies a favorable gain in free 
energy when bringing a hard sphere from the bulk to the wall. 



leading, after integrating over a half-space, to the contact 
relation, 



P 



2e' 



(23) 



where P w is the local pressure at contact with the wall. 
The notion of a local pressure can be somewhat mislead- 
ing, since the pressure is normally associated with bulk 
properties. It is only sensible to talk about a local pres- 
sure within the framework of a local density approxima- 
tion. 

Going back to the exact contact value theorem in 
Eq. (|14[) shows that the correct thermodynamic result 
is obtained only if P^ = ksTpb- Paradoxically, if one 
tries to improve on this result and go beyond the PB de- 
scription by including an excess free energy due to elec- 
trostatic or steric correlations using a LDA formalism, 
one obtains less accurate contact value density than if 
the correlations are neglected altogether. The unphys- 
ical contact value relation for the LDA treatment is a 
consequence of an unphysical density profile. The LDA 
approximation fails to capture the oscillating discrete 
structure of a fluid, see Fig. As seen in that fig- 

ure, the LDA approximations account for overcrowding 
through "density saturation", while a realistic descrip- 
tion shows overcrowding as layering of a fluid near a wall. 
It is tempting to assume that the "saturation" observed 
in a density for the LDA models represents a weighted 




FIG. 3. The density profiles near a hard charged wall obtained 
using the DFT and the two LDA methods: the MPB and 
the SPT. The parameters are: c s = 0.1M, a c = 0.4C/m 2 , 
As = 0.72nm, and a — 0.8nm. The system represents a 
symmetric 1 : 1 salt. 



rather than a physical density profile, 

dx'p{x')u{x' -x), (24) 



„LDA 



P 



where lo(x' — x) is some weight function, so that oscil- 
lations are completely smoothed out. Such conjecture, 
however, entails that tf> obtained from the LDA is also 
not physical but only a weighted quantity. We, there- 
fore, consider that all quantities produced by the LDA 
are what they are, and that the saturation of a density 
profile is an artifact of the LDA treatment. To test accu- 
racy of the LDA one should rather consider other quanti- 
ties such as ip, which is of direct interest to most problems 
of electrostatics. 

Below we look into details of two LDA models. As 
mentioned earlier, the SPT results are highly accurate 
for homogeneous fluids, producing the exact second virial 
coefficient. Furthermore, the third and higher order virial 
coefficients are very accurate, 



P P «* - (1 - 7? )3 - P 



1 + 4?7 



(25) 



On the other hand, the lattice-gas model (from which 
the MPB is derived) already fails at the second virial 
coefficient, 
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ig 



iog(i - n) = p 



l 



i 



(26) 



and is unfit for treating homogeneous fluids. 

On the other hand, when the LDA is applied to inho- 
mogencous fluids, it is hard to predict what will happen. 
In Fig. (j4]) we plot the contact potential for various mod- 
els. First we note that the difference between the PB 
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and the DFT results. The hnite size effects lead to an in- 
crease of the surface potential, and the larger the surface 
charge, the larger the discrepancy, since more counterions 
are crowded into the double-layer. Next we discuss the 
LDA models. The MPB model underestimates while the 
LDA based on the SPT overestimates the DFT result. 




CO. 



a c [C/m z ] 



FIG. 4. The contact potential as a function of the surface 
charge. The parameters are: c s = 0.1M, As = 0.72nm, and 
a = 0.8nm. The system represents a symmetric 1 : 1 salt. 

We next consider the quantity A w as defined in 
Eq. ([TT|). As shown in Fig. © the MPB reproduces this 
quantity better than the STP. The advantage of the MPB 
model is its analytical tractability, and according to this 
model A„, yields a very simple expression, 



2^-3 



I2e 



(27) 



with A w parabolic in a c . 




a c [C/m 2 ] 



FIG. 5. The reversible work of bringing an uncharged test 
particle from the bulk to the wall as a function of the surface 
charge. The parameters are: c s = 0.1M, As = 0.72nm, and 
a = 0.8nm. The system represents a symmetric 1 : 1 salt. 



the contact value theorem is 
16 

fivP + — X* B a* 2 = Tj b cosh t 



(28) 



■KG' 



/6 



where <ft = e/3ip is the reduced potential, and v 
is the volume of a single fluid particle. 

According to the LDA, A w is a function of two inde- 
pendent parameters, r\b and X B a* 2 . Later, as we inves- 
tigate A w for various parameters using the DFT, we do 
not find this reduced description to be confirmed, and we 
conclude that it is an artifact of the LDA. Furthermore, 
this implies that the scaling of the contact potential ij) w 
according to the LDA description is described by two 
independent parameters. The two parameter scaling of 
the LDA can be ascertained from the contact relation in 
Eq. (E2D. 



IV. SCALING OF THE QUANTITY A w 

All the calculations in this work are carried out for the 
1 : 1 symmetric salt, so that the bulk number density 
is pb = 2c s , where c s denotes the bulk salt concentra- 
tion. In this part of the study we investigate the scal- 
ing of the quantity A w . A w is a function of three di- 
mcnsionless variables: rjf, — ira 3 pb/6, cr* = no c o 2 / (4e 2 ), 
and X B = X B /a, corresponding to volume, surface, and 
length, respectively, a* represents the packing fraction of 
counterions projected on a 2d plane. The 2d solid phase 
transition is at r/2d > 0.72, and the maximum possible 
packing fraction is at X]2d = 0.90. If a* < 0.90, potentially 
all the counterions can collapse onto the charged surface, 
if attraction is sufficiently strong (although a* > 0.72 
would involve a phase transition as well) . If, on the other 
hand, a* > 0.9, collapse onto a single plane is prevented 
and a second fluid layer must form. In the reduced units, 



A. dilute limit 

To simplify things we first consider the dilute limit, 
so as to suppress the contributions from the bulk elim- 
inating the dependence on r/b, lim^^o A w = A^ = 
A^(cr*, X* B ). The main results are in Fig. © and Fig. 0. 
The figures show the scaling behavior of A w . The general 
observation is that the scaling with the surface charge is 
stronger than that with the Bjcrrum length. This can be 
understood since the surface charge controls both the in- 
teractions between the surface and counter charges and 
the counterion concentration. On the other hand, the 
Bjerrum length regulates only the strength of the ionic 
interactions. 

Our aim is to propose a fit which reproduces the 
scaling of A w . The desirable procedure would be to 
suggest an expression based on a physically motivated 
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The accuracy of this fit is demonstrated in Figs. (|6]) and 
Fig. (|7|) in comparison to the DFT data points. We also 
compare this formula with the simple expression from the 
MPB model in Eq. (f2~T)) which in reduced units reads 



A 1( V 

= y c 



: \ * 



(30) 



There is a curious similarity between the square term in 
our fit = 5.53 and the coefficient 16/3 in the MPB ex- 
pression. If this similarity is more than coincidental, it 
would suggest that although reduced, the MPB descrip- 
tion accurately captures some aspects of the hard-sphere 
interactions. 



B. finite concentrations 



FIG. 6. A„ as a function of the surface charge. The packing 
fraction is fixed at rjb = 0.0040 (the dilute limit). The points 
are the DFT results and the lines are given by Eq. (|29[) . 
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FIG. 7. A w as a function of the Bjerrum length. The packing 
fraction is fixed at 771, = 0.0040 and represents (the dilute 
limit). The points are the DFT results and the lines are given 
by Eq. ([29]). 



We now turn to investigate the contributions of finite 
concentrations. We suggest the following formula, 



Ai /f — A.., 



log 



P 



k B Tp b 



= A° w -6(r, b ,X B ), (31) 



so that A w as a function of a* is merely shifted by the 
quantity 6. This fit is suggested by Fig. ([5]) where we 
compare A w for finite rjb with that for a dilute limit. In 
the relevant regime A w as a function of a* is offset by a 
value independent of a*. We find the following accurate 
fit to the function 9, 



1.7577b- 



% 2 



10.8A1 



1 



0.21 + X* f 



(32) 



We subtract — \og{P /ksT pt,) from A w so that the quan- 
tity A w is always positive, see Eq. (flU)) and Fig. @. To 
ensure that < A° , so that A w > for all values of rfc, 
we apply the form, 



A?„e-^. 



APPLICATIONS 



(33) 



model. Unfortunately, an accurate physically motivated 
fit could only be found for a limited range. To ac- 
curately parametrize the entire range we had to resort 
to a fitting procedure where the choice of a functional 
form was based on simplicity and convenience. The 
following functional form was found convenient, A^, = 
a(a*){X* B /[l + b(a*)\* B ]} a . Fully parametrized, this form 
reads 



A° = (2.28ct c * + 5.53ct* 2 ) 



1 + (1 - 1.210-* + 0.375cr* 2 ) X* B 



3/4 



(29) 



We now have the complete expression for A w , 



A w = -log( ) +A W , 



(34) 



where A w is given in Eq. (f3"3"|) . To obtain the contact 
potential, we insert A w into Eq. ([28)) to get 



(35) 



cosh^= e Ml 1 16 AW * 2 



3 v(3P 



In Fig. © and Fig. (fit)]) we compare the exact DFT 
results for the contact potential with the ones obtained 
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FIG. 8. A w as a function of the surface charge for different 
Bjerrum lengths and packing fractions. 
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FIG. 10. Surface potential as a function of the particle diam- 
eter. The parameters are: a c = 0.2C/m 2 and As = 0.72nm. 
The solid lines correspond to our analytical expression and the 
dashed lines to the MPB model. Since the hard sphere fluid 
freezes at rjb — 0.49, we interrupt the curve for c 3 = 0.5M. 



The standard PB result gives 



Cqc = eKCOsh 



which is the, so-called, Gouy-Chapman capacitance, k = 
\/8ttXbCs is the inverse Debye length. At low (j) w this 
formula works well, but as <p w increases, steric repulsion 
lowers the capacitance as countcrions arc prevented from 
accumulating near the electrode, leading to a "camel- 
shaped" curve of Fig. (fTTj) . The scaling function derived 
in the present paper allows us to easily calculate the dif- 
ferential capacitances, which once again are in a very 
good agreement with the DFT. 



FIG. 9. Surface potential as a function of the surface charge. 
The parameters are: As = 0.72nm and a = 0.8nm. The 
solid lines correspond to the analytical expression based on 
Eq. (J35J) and the dashed lines to the MPB model. 

using the analytical scaling function Eq. (|35[) and with 
the predictions of the MPB model. The scaling form 
derived in the present paper agrees very well with the 
contact potential obtained using the exact DFT, even at 
densities all the way to the solid phase transition. For 
largest particle sizes, Fig. (fTOj) shows that the scaling 
function slightly disagrees with the DFT results. This 
happens when a* > 2 and a third layer of condensed 
counterions starts to form. 

A quantity that is of particular interest for electro- 
chemistry is the differential capacitance of the double- 
layer [Hf, 

C = — 

dip w 



VI. CONCLUSION 

We have investigated the effects of steric interactions 
on the ionic distribution and electrostatic properties near 
a charged hard wall. We have specifically focused on 
the reversible work required to bring a hard spherical 
particle (with no electric charge) from the bulk to the 
wall, A^. The quantity A w is nonlocal and it controls 
the steric effects in the double-layer. Typically A w > 0. 
This implies that the overcrowding near a wall causes the 
depletion of counterions away from the wall, leading to 
the increased contact potential. The less common case, 
A w < 0, represents the situation where counterions are 
forced from the dense bulk toward the wall. This effect 
decreases the value of the contact potential. 

We have also studied the scaling behavior of A w calcu- 
lated using the nonlocal DFT theory. We expected that 
the careful look into the scaling behavior can reveal a 
simple underlying mechanism. Unfortunately, no simple 
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FIG. 11. Capacitance versus the surface potential. The pa- 
rameters are: As = 0.72nm and a = 0.8nm. The solid lines 
correspond to the analytical expression, the dashed lines are 
the result of the MPB model, the symbols are the DFT data 
points, and the dotted lines are the PB results, Cgc- The 
agreement between the DFT and the simulation for C was 
also confirmed in [25|. 

model could cover the entire domain. Instead we pro- 
posed a functional form to fit the DFT data points for 
A w . Such a fit allows for simple and accurate calculation 
of electrostatic properties at contact with the wall with- 



out resorting to numerical schemes of an entire region. 

Finally, we have explored the discrepancies between 
the LDA and exact description. The LDA description 
shows that A w depends on two free parameters, in com- 
parison to the three parameters given by full physical 
description (DFT). Another artifact of the LDA is the 
" saturation" seen in a density profile rather than oscilla- 
tions due to layering. This saturation effect leads to un- 
physical contact value relation. For various LDA models, 
we find that the MPB based on the lattice-gas model is 
more accurate than the LDA based on the SPT theory, 
which for treating a uniform fluid is extremely accurate. 
This indicates that some cancellation of errors is at work. 

We remind that our analysis is limited to the weak- 
coupling limit, where electrostatic interactions are ac- 
curately captured by the mean- field treatment. To ac- 
curately account for the correlations together with the 
hard-sphere interactions remains a great challenge. It is 
our hope to treat the strong-coupling limit of the hard- 
sphere ions in the future. 
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